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Abstract 

We investigate the manifold Ai of quadratic forms in n variables having a 
multiple eigenvalue. We present an algorithm to derive a minimal-degree equation 
system for Ai, and give its result equations for n = 3. We prove that 1) At is 
prime and has co-dimension 2, 2) each matrix in Ai having n — 1 of different 
eigenvalues is a regular point on the surface At, 3) in the case n = 3, the set 
of singular points on Ai is the set of scalar matrices. We give a geometric 
description of Ai in a neighborhood of each regular point: a fibration over a plane 
with the fiber being an orbit by conjugations by SO(n). For n = 3, Ai is also 
described as the straight cylinder over Aio, where .Mo is the cone over the orbit 
of the diagonal matrix diag(l, 1, —2). We analyse certain interesting metrical and 
algebraic properties of this orbit, which also occurs a diffeomorphic image of the 
projective plane. These results simplify, generalize and complete the results given 
in some previous works on this subject. 

Keywords: quadratic form, orthogonal change, discriminant, dimension of manifold, 
primality of manifold, singular point. 



1 Introduction 

About the word "fibration" 

We use it in the two meanings: 1) fibration over a plane considered modulo isometry, 
2) fibration over a plane considered modulo diffeomorphism. 

Let us denote: 

Sym(n) = Sym(n, M.) is the space of quadratic forms — that is symmetric matrices 
(of size n) over the field 1R of real numbers (dimSym(n) = n(n + l)/2), 

charPol(X)(A) is the characteristic polynomial of a matrix X of Sym(n), 
Ai = {X G Sym(n)| discr X = 0} is the surface in Sym(n) defined by discriminant 
of charPol X. In other words: Ai is the variety (manifold, surface) of all matrices in 
Sym(n) that have any multiple eigenvalue. 
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The problem is: to describe the geometry of the algebraic variety Ai. 

About generality: our conclusions are valid for each dimension n. But for n < 3, 
Ai is the line of scalar matrices, and all the below constructs occur trivial for this case. 
Therefore we consider below only the case of n > 3. 

It is known that each matrix in Sym(n, M) has n real eigenvalues, some of which may 
coincide. 

We call the multi-set of eigenvalues the set of eigenvalues for a matrix in Sym(n) 
— together with their multiplicities. 

The surface Ai is given by a single equation discr(X) = 0, where the polynomial 
discr(X) has the matrix elements Xij as variables, has degree 2n and many monomials 
- even for n = 3. A knowledge about the surface Ai has an important application, 
especially for the case of n — 3. Therefore there were done several investigations for 
finding various equation systems for Ai. Thus, for n = 3, it occurs (|Illj. 1992) that 
discr(X) is the sum of squares of several simpler polynomials. For example, it is visible 
from this that there may be needed extra considerations to prove that co-dimension of 
Ai is less than 3. 

It is probable that a large part of this problem has been solved long ago by classic 
scientists, especially for n = 3. But it is unlucky that I failed to find sufficient references. 
Currently, I see the following literature. 

In |Ar] of 1999 (I have seen only Russian version) it is written in the Section 6 about 
relation of this problem on quadratic forms to the von-Neumann-Wigner theorem about 
intersection of electronic levels. It is also written there (I translate from Russian): 
'As we have seen, the theories of eigenvalues of quadratic, Hermitian, and hyper- 
Hermitian forms are actually the theories of hyperbolic cones in the spaces of quadratic 
forms" . 

Symmetric matrices in our problem are actually quadratic forms. 

|Ar] writes further, after Corollary 3, about co-dimension d + 1 for "the variety of 
forms with characteristic roots of non-minimal multiplicity". Our case is of the real 
field, so, we could understand the above statement, — together with the statement in 
the next paragraph, — so: "co-dimension of Ai is 2". 

But it is difficult for me to compile from this paper a simple proof for this statement. 
We provide here a direct proof for this, together with the solution for other parts of the 
problem. 

In the papers [ffi EI] (of 2004 and 2005) 

1) it is written that Wigner and von Neumann considered the variety Ai (for n = 3) 
and have provided certain in-formal reasons of why its co-dimension should be two, 

2) it is written that in [IkJ there is proved (for n = 3) that this co-dimension is not less 
than two, and this proof uses the result of the paper |I11] about decomposing discriminant 
into a sum of squares 

(we add that there also exists a trivial proof for this inequality), 

3) there are given (in jlk | ID:Ij ) the reasons of why in the case of n = 3 the dimension of 
Ai is 4, and it is given a certain discourse about the irreducible components. 

And these considerations include various transformations with explicit equation sys- 
tems for Ai. 

It this our paper we also deal with the case of arbitrary n, and apply a different 



2 



approach: the classical method with a linear Lie group of symmetries. 

In brief, our approach is as follows. All symmetric matrices are produced by 
conjugating diagonal matrices with operators from SO(n). A displace in the plane of 
diagonal matrices is orthogonal to a displace along this conjugation orbit. This provides 
a smooth parameterization for A4, except certain particular points. And it remains to 
describe the orbit for a diagonal matrix, with considering the two cases for the number 
of its different eigenvalues: n — 1 or less than n — 1. 

As to finding for Ai of a minimal (in total degree) system, this problem is solved by 
a simple and generic tool: combining of the Grobner algorithm for a polynomial ideal 
with diagonalization by the action of the 0(n) group. 

The results of this paper are as follows. 

• We simplify the conclusions and discourse. 

• We prove that Ai is prime and has co-dimension 2 
(Sections 121 l5l l5T4j) 

(the decomposition to primes suggested in [D:Ij is erroneous). 

• We prove that the matrices in Ai having the number n — 1 of different eigenvalues 
are regular points on Ai (Section 15.41) . 

For n = 3 we prove that scalar matrices and only them are singular points on 
M (Section EU). 

• We describe the global structure of the variety Ai in the n(n + l)/2 - dimen- 
sional space. For a neighborhood of a matrix of maximal spectrum, we describe a 
diffeomorphic parameterization as a fibration over a plane. 

For n = 3, we also provide a more definite description for Ai: the straight 
cylinder over Aio, where Aio is the cone over a diffeomorphic image of projective 
plane (Sections IU 15. 6p . 

• We describe certain metrical and algebraic properties of the orbit — as of a surface 
in Sym(3) (Section l5.2j) . 

• We give a simple algorithm for deriving a minimal (in total degree) equation system 
for Ai, — and also for any separate orbit, — and provide its output for the case 
of n = 3 (Sections l3l 15. ip . 

So far, we leave un-solved the following problems. 

1) Find whether it is true for n > 3 that each matrix in Ai having narrowed spectrum 
is a singular point. 

2) For n = 4, find minimal-degree equations for Ai (the algorithm given below is 
rather expensive in computation for n > 3). 

3) The structure of the surface Ai depends mainly on the structure of an orbit 
0(A, n) for a matrix A of maximal spectrum in A4. Similarly as for n = 3, it would 
have sense to find for n = 4, 5 some global description for 0(A, n), more definite than 
the one that we give in the sequel. 

To reduce the main text volume, some of not so interesting proofs in this paper are 
moved to Application, and the main part has a few references to Application. 
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2 Several definitions and preliminary constructs 

1. Dg denotes the space of diagonal matrices (it has dimension n), 

Seal denotes the line of scalar matrices. M is a conic set containing Seal. 

2. For a diagonal matrix D, denote D(i) = D(i,i). 

diagf ) denotes the diagonal matrix D having D{1) — ai, . . . , D(n) = a n . 

A diagonal-ordered matrix is a diagonal matrix D in which -D(l) < . . . < D(n). 

3. 0(n) is (an algebraic and linear) Lie group of orthogonal operators in M. n . 
SO{n) is (irreducible and connected) Lie subgroup of rotations in 0{n) (determinant 
= 1), it has dimension dimSO(n) = n(n — l)/2. 

4. X — > g ■ X ■ g^ 1 is the action of the coordinate (basis) change by an orthogonal 
operator g in a matrix (operator) X. 

We shall shortly call this conjugation action "change by g" , "changes by operators 
from SO(3)" , "changes from SO(3)" , and such. 

And for an orthogonal operator g, there holds g' 1 = g* (transposed matrix), and 
the expression g ■ X ■ g* represents the basis change in a symmetric bilinear form X. 

For g in O(n), we denote the action of the basis change as 

g c X= g c {X)= g-X-g*. 

Further: g~ c denotes (# _1 ) c , and l* = T(Z, 0) c denotes (for n = 3) the 
operator of the change by the rotation operator T(l, <fi) for the axis / and angle 0. 

In the case of n — 3, we call 1-orbit of a matrix M by the axis / the orbit of M 
under the changes by the rotation operators around I. 

5. Denote MV the (conic) set of diagonal matrices in M; 

for 1 < i < j < n denote IL j the plane in Dg defined by the equation D(i) = D(j). 
Evidently, AiT> is union of the number C\ of (n — 1) -dimensional planes Iljj. 
For example, for n — 3, II = LT12 is all the matrices of the kind diag(A, A,/i), while 
LT13 and 112,3 are obtained from II by permutational changes from 5*0(3). 

6. Below, the word "orbit" means: the orbit of a matrix in Sym(n) under the action 
of changes by the operators from SO(n); Orbit (A) denotes the orbit of a matrix A. 

7. Fix an orthonormal basis Bas = {ei, . . . , e n } in M. n , and let us consider the 
matrices in Sym(n) as representations in Bas of symmetric bilinear forms. We also 
speak of operators from SO(n) given by permutations on the set Bas, or, for example 
(in the case of n — 3), being rotations around the axis e*. 

8. trace(X) for a matrix X is the sum of the elements on its main diagonal. 

Generally, our simple discourse bases on 1) known facts of the linear Lie group 
theory, 2) classical 
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Theorem (DO): for any real symmetric bilinear form there exists an orthogonal 
change which brings this form (matrix) to a diagonal matrix 
(see, for example, |VWj . paragraph 90). 

Lemma. For such a diagonalization there are sufficient the operators from SO(n). 
Indeed, if g c diagonalizes X, and detg = — 1, then the composition g2 of g with the 
permutational operator for (ei,e2) belongs to SO(n), and g\ diagonalizes X. 

Several more definitions 

1) The stabilizer St (A) of a matrix A from Sym(n) is the set of operators g, which 
changes preserve A. It is a smooth subgroup. 

2) The width of the spectrum for a matrix from Sym(n) is the number of its 
different eigenvalues. 

For a set X in Sym(n), any matrix from X having maximal in X spectrum width we 

call a matrix of maximal spectrum (in X), 

and all the rest are called the matrices of narrowed spectrum. 

3) Denote AiH the set of matrices of maximal spectrum in Ai — that is having 
the spectrum width n — 1. This is a topologically open set in Ai. 

4) Denote EL^ = IL^ n AAU. 

The following; three statements are evident. 

(1) The spectrum width of a matrix in Ai is the spectrum width of its diagonal form. 

(2) For a matrix in Sym(n) having all different eigenvalues, the orbit dimension is 
dimSO(n), 

matrices of maximal spectrum in Ai have orbits of maximal (for Ai) dimension, equal 
to dimSO(n) — 1. 

(3) If a matrix in Sym(n) has the spectrum width n — 1, then its stabilizer is a single- 
parameter subgroup of rotations of the plane L(u, v) — for each two non-collinear vectors 
u and v belonging to the same multiple eigenvalue. 

Example: for n = 3, the diagonal matrices in AiH. are exactly the matrices of 
kind diag(A, A,/i) with A 7^ fi, and also the two families obtained from this one by 
permutations on the main diagonal. 

s-metric 

Define an Euclidean metric on the space Sym(n): sum of squares of the elements of a 
matrix X. We denote this quadratic form sQuad(X) 
and call it s-metric. 

In this paper, the expressions "distance", "angle", "orthogonality", "isometry", 
"orthogonal projection", "circumference", "sphere", "circumference center", "radius", 
"surface diameter", "straight cylinder" applied to points and subsets in Sym(n) 
are understood in the sense of s-metric. 

Lemma S (probably, known). 

For any matrix X from L(n) and any orthogonal operator g, it holds the equation 
sQuad(X) = sQuad(g • X) = sQuad(g c X). 

That is: the changes from 0(n) preserve s-distance and s- angles in L(n). 
Proof. An orthogonal operator preserves s-square for each column-vector in X. Hence, 
sQuad(gX) = sQuad(X). Similarly, the right-hand side multiplication by g* preserves 
s-square of each row. Therefore sQuad(g c X) = sQuad(X). 
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Lemma 1. 

(1) Any change from SO(n) maps M. on itself. 

(2) Scalar matrices are the only fixed points of the action of SO(n) on Sym(n). 

(3) M. is the union of mutually non-intersecting orbits. 

(4) Each diagonal matrix in M. is brought by some change from SO(n) to an ordered 
diagonal, (that is all the results for permutations applied to the main diagonal of a 
matrix in AiT> are represented by some changes from SO(n)). 

(5) Two matrices in M belong to the same orbit if and only if they have the same 
eigenvalue multi-set. 

(6) For each matrix in Ai, its orbit has exactly one diagonal-ordered matrix. 

(7) 7.1. For each real number t and symmetric matrix M, it holds 
Orbit(t-M) = t ■ Orbit (M). 

7.2. For n = 3, orbits of any two matrices in MV, differ from each other in a shift 
at some scalar matrix and a homothety by some non-zero factor 
(and we use that this transformation preserves angles). 

(8) 8.1. Discriminant of a symmetric matrix does not change with adding to this 
matrix of any scalar matrix. Any shift by a scalar matrix maps the surface Ai on itself. 

8.2. For each real number s, intersection Ai s of Ai with the hyper-plane trace(X) = 
s is a surface mapped on itself by the changes from SO{n). 

Ai is a straight cylinder over Mq having the line Seal as element. 

Proof. The statements (1), (2) and (3) are known and evident. 

Let us prove (4). Let the searched permutation on the main diagonal is presented 
by an orthogonal operator g. D has at least two equal elements on the main diagonal. 
Hence, the permutational operator g 2 for the corresponding basis vector pair has deter- 
minant —1. Therefore, either g belongs to SO{n) or the composition g • g 2 belongs to 
SO (n) and its change brings D to a diagonal-ordered form. 

Let us prove the statements (5) and (6). The change action preserves the coefficients 
of the characteristic polynomial. Hence, it remains only to prove the second part of the 
statement (5). Let A and B have the same eigenvalue multi-set. By Theorem DO 
and its additional lemma, A and B are diagonalized by some changes from SO{n). By 
the statement (4), they are brought further to a diagonal-ordered matrices. As these 
diagonal-ordered matrices represent the same multi-set, they are equal. 

This also proves the statement (6). 

Let us prove the statement (7). The statement (7.1) is evident. As to (7.2), it is, 
evidently, sufficient to prove it for a pair of diagonal matrices in AVH. And hence, it is 
sufficient to prove it for a pair of matrices (D,D 2 ), where D = diag(0,0, 1), and D 2 
is any matrix of kind diag(A, A, jj) with A^/j. 

D 2 is obtained from D by adding of some scalar matrix Sc = diag(a, a, a) and 
by homothety by some non-zero number factor b: (D + Sc) ■ b = D 2 . For such 
representation, it is sufficient to put a = A/(/i — A), b = /i — A. Each change operator 
g c is an isomorphism on the algebra of square matrices. In particular, for each square 
matrix M there hold the equations 

g c (M + Sc) = (g c M) + (g c Sc) = (g c M) + Sc, g c (b ■ M) = b ■ g c M. 
Therefore with adding of a scalar matrix Sc to any symmetric matrix, the orbit shifts 
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at the vector Sc, and with multiplying M by a coefficient b, the orbit is homothetically 
multiplied by b. 

Let us prove the statement (8). 8.1: Adding a scalar operator \i-E to a symmetric 
operator M shifts the spectrum of M at the number /i. This preserves discriminant, 
because the discriminant depends only on the differences of eigenvalues. 

8.2: For each real number s, denote IT* the hyper-plane in Sym(n), defined by the 
equation trace(X) = s. Denote Ai s = Aid IT*. Any basis change operator preserves 
the matrix trace. Therefore each such plane is mapped on itself by the changes from 
SO(n) — as well as each restriction Ai s . 

Further, the surface Ai is union of the restrictions Ai s for all s. Adding of any scalar 
matrix s ■ E shifts the surface Ai along itself — because (by the statement (8.1)) this 
adding preserves discriminant. This shift maps the plane ITg on the plane II*, and maps 
the restriction Ai on Ai s . Hence, Ai is the cylinder over Aio having the line Seal as 
element. The line Seal is orthogonal to the plane ITg, as one can see from the s-product 
of a scalar matrix by a matrix of zero trace. Hence, this cylinder is straight. 

The lemma is proved. 

Lemma AO (probably, known). 

Commutator of a diagonal and an anti-symmetric square matrices (of the same size) 
is a symmetric matrix having zero diagonal. 

Lemma Al (known). The counter-image Ai' of an algebraic set Ai for a polynomial 
map F : R n — > R m is an algebraic set. 

Lemma A2 (known). For a surjective polynomial map from an algebraic set Ai 
onto an algebraic set Ai', if Ai is prime, then AM is prime. 

Simple proofs for the Lemmata AO, Al, A2 are given in Application. 

Theorem "Primality". 

The discriminant surface Ai is algebraically prime 

(that is it is not a union of any two algebraic sets which do not contain each other). 
Proof. Consider the map DO : II x SO(ti) -> M, 

DO((X, Hi, . . . , Hn-2),g) = 9 C diag(A, A, . . . , /i n _ 2 ) 
(no restriction for the numbers in the first argument). This is a surjective polynomial 
map of algebraic surfaces, where the support domain has one dimension more than the 
image Ai. The algebraic surface n x SO(n) is prime. Hence, by Lemma A2, Ai is 
prime. 

To be fully rigorous, we need to ensure that the image of DO is Ai. By Theorem 
(DO) and its additional lemma, replacing n with AiT> makes this map surjective. By 
Lemma 1 (4), all AiT> is obtained from n by changes from SO(n). 

The theorem is proved. 
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3 General algorithm for minimal— degree equations 
for M 



Initially, Ai is denned by a single homogeneous form discr of degree 2n. It is natural 
to search for possibly smaller-degree equation system for this surface. Let us formulate 
a generic algorithm for this. 

The surface Ai is defined by the ideal I{Ai) of all polynomials in the variables Xs 
which are zero on Ai. Let us call this ideal a complete system for Ai. 

Fix the grading by the total degree on this polynomial algebra. There is known the 
algorithm of the Grobner basis |Buj . Given any finite basis for an ideal in a polynomial 
ring (over an Euclidean coefficient ring), this algorithm produces a finite basis for this 
ideal minimal by the given grading, and, in a certain sense, canonical. In our case, denote 
this Grobner basis gs. The property of the method is so that if the ideal I {At) contains 
any polynomial of total degree less than d, then gs contains a polynomial of total degree 
less than d. Therefore, the problem of defining of the surface Ai by equations of possibly 
small degree has its generic and algorithmic solution — if only there is found any finite 
basis for a complete system for the surface. For example, {discr} is not a complete basis 
in our case, as we shall see further. 

In general, our method is as follows. If a map F : R m — > Ai is surjective, then 
this defines the maps X(i,j) from M m for each position in the matrix X. Then, 
the ideal Rel of all the algebraic relations between the maps {X(i,j)\l < i < j < n} 
is just a complete equation system for Ai. It remains to formulate the case when a 
parameterization F (not necessary injective, or smooth) allows an algorithm for finding 
of a finite basis for Rel. 

In our case, we provide as F a certain polynomial parameterization. Namely, build 
a diagonal matrix 

D = diag(A,A,/ii,...,/i„_ 2 ) 

and consider its diagonal elements as variables for polynomials with coefficients in M. 
Denote Eigs the set of these variables. Further, let the matrix Y of size n x n consist 
of the elements 

Ys = {y id | l<i,j< n}. 

We denote the variable set: Zs = YsUEigs, and we shall say that the variable set 
consists of n 2 of o-variables and n — 1 of eg-variables. Consider the matrix product 
X = Y ■ D ■ Y*, where the factors are considered as matrices over R[Zs]. 

Evidently, a generic diagonal matrix of Ai is presented by the above D and also 
by the matrices made from D by permutations on the main diagonal. By Theorem 
DO, each matrix in Ai is diagonalized by some orthogonal change. This gives us the 
following generic matrix for Ai: 

X = Y-D-Y*. 

Here the elements of these matrices are considered as elements of the quotient-ring 
QZ of the algebra R[Zs] by the orthogonality conditions for the operator Y. These 
conditions form the set 

OrtEs = {row(z, Y). row(j, Y) = 5(i,j) | 1 < i < j <= n} 
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-number n(n + l)/2 of quadratic equations. Here the expression row(i, Y) denotes 
taking the row No i from a matrix, and the symbol "." denotes scalar product of two 
rows. 

The above generic matrix X uses only one matrix from the set of diagonal-permutati- 
on results for D, because all the diagonal permutation results in D are expressed by 
changes from SO(n). 

The generic matrix X (over QZ) is symmetric. The elements X(i,j) (conjugation 
classes) have representatives as polynomials in M[Zs]. Now we see that the ideal Rel of 
algebraic relations for X(i,j) is a complete equation system for Ai (this ideal consists 
of polynomials in M[Xs] in the variables Xs = 1 < i < j < n}). 

Further, there is known a simple algorithm for finding of a canonical basis Rels for 
the ideal Rel. It is described in |GTZj (Corollary 3.2), and it bases, in its turn, on the 
Grobner basis method. This algorithm produces such a basis for the ideal of algebraic 
relations that is minimal in the grading and canonical. For example: if there exists a 
polynomial of total degree less than three which is zero on the surface Ai, then the 
returned basis Rels also must have such a polynomial. 

About the coefficient domain The problem formulation mentions the surface in a 
real space and its equations over the real number field. And even for a matrix in Ai 
having rational elements, its eigenvalues are, generally, algebraic numbers. 

Nevertheless: discriminant in our problem has integer coefficients, for natural rea- 
sons, the obtained system for Ai has integer coefficients, and all the intermediate com- 
putations need only to operate with polynomials having rational coefficients. 

About computer calculation 

The algorithm presented above is important independently on a computer device. 
For example, it can be evaluated by hand for each given value of n (though, it is difficult 
to compute for n > 3). On the other hand, computer evaluation is also interesting and 
helpful. We applied for this our common system for computer algebra called DoCon. 

Case n = 3 

For this argument, we have obtained the result of the above algorithm. It also can 
(with some effort) be evaluated by hand, though we skip this exercise. The result is the 
system Rels for the surface Ai containing seven homogeneous cubic equations and one 
more equation of degree four. 

Note: Rels does not contain any polynomial of total degree less than three. Hence, 
the surface Ai is essentially non-quadratic. This is due to the following property of a 
Grobner basis: if a non-zero polynomial is zero when restricted to Ai (belongs to I(Ai)), 
then its leading monomial is a multiple of leading monomial of some polynomial in the 
basis Rels. 

Further, if Rels has any polynomial that has some power belonging to the ideal of 
the rest of the system, then this polynomial can be removed. And let us repeat this 
simplification while such a polynomial is found in the current system. There exists a 
certain simple algorithm which detects the above relation between a polynomial and 
a finite set of polynomials. This process produces a reduced system defining the same 
surface. In our example, it shows, again, four homogeneous forms (it prints the left-hand 
sides of equations, while the right-hand sides are zero): 
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RelsS = 



{xl2~2*x23 -xl2*xl3*x22 +xl2*xl3*x33 -xl3~2*x23, 
xl2*x23*xll -xl2*x23*x22 -xl3~3 +xl3*x23~2 -xl3*xll*x22 +xl3*xll*x33 

+xl3*x22~2 -xl3*x22*x33, 
xl2*xl3*xll -xl2*xl3*x22 -xl3~2*x23 +x23~3 -x23*xll~2 +x23*xll*x22 

+x23*xll*x33 -x23*x22*x33, 
xl2"2*xll -xl2~2*x22 -xl3"2*xll +xl3~2*x33 +x23~2*x22 -x23~2*x33 -xll~2*x22 

+xll~2*x33 +xll*x22~2 -xll*x33~2 -x22"2*x33 +x22*x33~2, 
xl2~3 -xl2*x23~2 -xl2*xll*x22 +xl2*xll*x33 +xl2*x22*x33 -xl2*x33~2 

-xl3*x23*xll +xl3*x23*x33> 

Equations of this kind have been derived earlier by N. V. Ilyushechkin 

Certain four cubic equations IlEs(3) for Ai have been obtained manually and proved 
by N. V. Ilyushechkin for the case of n = 3, and we hope that the paper about this will 
be published in 2010 or in 2011. 

Our equation system RelsS is obtained by computer evaluation, after which we have 
compared these two systems and verified that the sets IlEs(3) and RelsS are equivalent: 
square of each element of one system belongs to the ideal of another system. 

Conclusion We 

1) present a generic algorithm for minimal-degree equations for Ai, 

2) give its computer result (Rels, RelsS) for n = 3, 

3) use the Ilyushechkin's equations to verify our result for n = 3. 

Equations for Mo 

By Lemma 1, M is a straight cylinder over Mo. Hence, it is sufficient to describe 

the surface Aio - the restriction to the plane xll + x22 + x33 = 0. xll is 

eliminated, and this results into the system 
MOeqs = 

{xl2*x23*x22 +(1/2)*x12*x23*x33 +(l/2)*xl3~3 -(1/2)*x13*x23~2 -xl3*x22~2 

+(1/2)*x13*x22*x33 +(1/2)*x13*x33~2, 
xl2~2*x22 +(1/2)*x12~2*x33 - (1/2) *xl3~2*x22 -xl3~2*x33 - (1/2) *x23~2*x22 

+(1/2)*x23~2*x33 +x22~3 +(3/2) *x22~2*x33 -(3/2)*x22*x33~2 -x33~3, 
xl2~2*x23 +(3/2)*xl2*xl3*x33 - (1/2) *xl3~2*x23 -(l/2)*x23~3 +x23*x22~2 

+(5/2)*x23*x22*x33 +x23*x33~2, 
xl2~3 -xl2*x23~2 +xl2*x22~2 +xl2*x22*x33 -2*xl2*x33~2 +xl3*x23*x22 +2*xl3*x23*x33} 

— again, four cubic forms, but the variable xll is eliminated. 
Checking computer evaluation 

How reliable is computer evaluation for this task? As a rule, large programs have 
errors, the same is with the very electronic schemes of computers. 

In our example, there are possible, at least, the following relatively easy checks. 

1. Discriminant must reduce to zero by the basis Rels (belong to /(Rels)). 

Each member of Rels has some degree which is a multiple of discriminant (we expect, 
square, or degree four is sufficient). 

2. Substituting into Rels of the values x it j for the matrix D = diag(A, A, /ii, . . . , \i n -'i) 
must produce zero polynomial in the variables A, fj,i, . . . , fi n -2 — this is easy to check. 

3. The matrix X (over the ring QZ) must be symmetric, 
its trace must be 2 A + fj,i + . . . + /i n -2, 

its determinant must be A 2 • ji\ ■ . . . ■ fi n -2- 
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4 Geometry of discriminant surface in Sym(3) 



In the case of n — 3, we have a more definite description of this surface. 

(DS3) of the structure of the discriminant surface M. in Sym(3, M). 

(1) M. is a prime algebraic variety of co-dimension 2. 
Its singular points are scalar matrices and only them. 

M. is defined by four cubic forms, and there does not exist a non-zero polynomial of 
total degree less than three which is zero on M.. 

M is union of restrictions M S) where the restriction for s is expressed by the equation 
trace(X) = s. 

M. is the straight cylinder over M.q having the line parallel to Seal as element. 
Therefore, it is sufficient to describe the surface M.q — and its description is as follows. 

(2) A three-dimensional surface M.q in a five- dimensional space has the following 
structure. 

2.1. M.q is the cone, with zero as vertex, over the orbit of the diagonal matrix 
diag(l, 1, —2); this two-dimensional orbit resides on the four-dimensional sphere having 
center in zero. 

2.2. The orbit of each non-scalar matrix in M. is a two-dimensional and algebraical- 
ly-quadratic surface residing on a four-dimensional sphere. It is diffeomorphic to the 
projective plane, and we provide a simple formula for this injection. 

The orbit orthogonally intersects the plane Dg of diagonal matrices, and it has exactly 

three common points with Dg. 

It is not contained in any four-dimensional plane. 

Also this orbit is defined by four quadratic equations. 

2.3. The diameter of the orbit of any non-scalar matrix in M having eigenvalues A 
and /i is |A — 

(3) The map DO : U x SO(3) ^ M, DO((X, //), g) = g c diag(A, A,//) 

(with no restrictions on A and //) is a surjective polynomial map of algebraic surfaces 
— from five-dimensional one onto four-dimensional. 

(4) A natural smooth parameterization for the surface .Mo\{0} can be presented 
as a direct product of the motion along the element line of this cone by smooth param- 
eterization for the orbit of the intersection point of this line with the sphere. 

The orbit of each matrix in Ai is obtained from the orbit of D — diag(l, 1, —2) by 
shift at some scalar matrix and by multiplying by some real coefficient. 

(5) Intersection of Mo with the plane of diagonal matrixes consists of the three 
different lines intersecting in zero. These lines have the following parameterization: 
diag(A,A, -2A), diag(A, -2 A, A), diag(-2A, A, A). 

Half of the proof for this theorem is given by the general theorem of the next Section. 
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5 Geometry of discriminant surface in Sym(n) 



Below, II denotes the plane IIi i2 of dimension n — 1 in the surface MT>. 

Theorem (DS) on the structure of discriminant surface Ai in Sym(n, R). 

(1) Ai is a prime algebraic variety of co-dimension 2. 

Ai is union of the restrictions Ai s , where the restriction for s is defined by the 
equation trace(X) = s. 

Ai is the straight cylinder over Aio having a line parallel to Seal as element. 
Matrices in Ai having n — 1 different eigenvalues are regular points on Ai. These 
matrices form an open and everywhere dense set AiH. in Ai. 

We give a simple algorithm (expensive to perform for n > 3) for finding for each given 
n of polynomial equations for Ai having minimal total degree. 

(2) Orbits for matrices in Ai are classified by the eigenvalue multi-set, or — by a 
unique diagonal-ordered matrix in the orbit. 

(3) AiH. has the following structure. 

3.1. AiH is a smooth fibration over the (n— l)-dimensional plane U h , with the fiber 
being an orbit — smooth compact surface of dimension dimSO(n) — 1. 

Each orbit in AiH orthogonally intersects the plane Dg of diagonal matrices and has 
exactly one common point with each sub-plane in Dg, conjugated with II by changes 
from SO(n) (there are number of C\ of such sub-planes). 
Orbit of a scalar matrix is a single point. 

3.2. Orbit of each matrix in AiV, is a smooth, compact, algebraic, connected surface of 
dimension dimSO(n) — 1 = (n(n — 1)/2) — 1. It resides at the sphere in the hyper-plane 
orthogonal to the line of scalar matrices. We give a simple algorithm (expensive to 
perform for n > 3) deriving from the given eigenvalues of a matrix A in AiH equations 
for the orbit of A having minimal total degree. 

3.3. Diameter dm of the orbit of a matrix in AiTi having the eigenvalue multi-set 
{Ai,...,A n }, is bounded as maxi — < dm < 2d(D,Scal), where d(D,Scal) 
is the distance from the matrix D = diag(Ai, . . . , A n ) to the line of scalar matrices. 

(4) The map DO : II x SO{n) -> M, 

DO((\, /ii, . . . , Hn-i),g) = g° diag(A, A, . . . , // n _ 2 ) 

(no restriction for the numbers in the first argument) is a surjective polynomial map 
of algebraic surfaces, where the support domain has one dimension more than the image 
M. 

(5) A smooth atlas on M.1-L can be naturally built as follows. 

From each point D = diag(A, A, fj,i, . . . , \i n -2) in AiTi there flows out the orbit Orbit(-D) 
of dimension dimSO(n) — 1. A smooth parameterization for this orbit in the neighbor- 
hood of D is given by "quotientation" of the action of SO(n). Namely, the stabilizer 
St(D) is a single-parameter subgroup of rotations of the bi-dimensional plane L(e±, e 2 ). 
This orbit is locally diffeomorphic to the homogeneous space SO' of conjugation classes 

by st(/;). 

A smooth parameterization for a neighborhood of each point of this orbit is by transla- 
tion to this neighborhood by some change from SO(n) of the parameterization near D. 
A total smooth parameterization for Ail-i is the product of the parameterization for the 
orbit of D and motion along a certain plane II in AiT>. 
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Now, let us proceed with the proof. 

First, recall that the statement (4) and the primality statement are proved earlier, 
in Theorem "Primality". 

5.1 Investigating an orbit by orthogonal changes 
Lemma 2. 

(1) If A is a non-scalar matrix in Sym(n), then its orbit intersects orthogonally the 
plane of diagonal matrices. 

(2) Orbit of any matrix in Sym(n) resides in the hyper-plane orthogonal to the line 
of scalar matrices. 

(3) Orbit of any matrix in AiH is a smooth, compact, algebraic, connected surface 
of dimension dimSO(n) — 1 = (n(n — l)/2) — 1. 

Dimension of an orbit for a matrix A in M. is for a scalar A, dimSO(n) — 1 - for 
a matrix of maximal spectrum, a number < d < dimSO(n) — 1 in other cases. 
Proof. 

(1) : Due to the symmetry by the group action and due to diagonalization, it is 
sufficient to prove this statement for an orbit of a diagonal matrix D. Also it is sufficient 
to prove the statement (A): each tangent vector for the orbit of D in D is some matrix 
M having zero diagonal. Then, by definition of the quadratic form sQuad, it will 
follow that sQuad(Di, M) = for each diagonal matrix D\. 

Now, prove (A). The tangent space T to the orbit of D in D consists of all the 
commutators of D with operators from the Lie algebra so(n), and this algebra consists of 
anti-symmetric matrices. Therefore, by Lemma AO, each matrix in T has zero diagonal. 
The statement (1) is proved. 

(2) : suppose the matrix A in Sym(n) is not scalar, and II* is the hyper-plane 
orthogonal to the line Seal and containing A. Let Aq be the base of this projection of A. 
By Lemma S, orthogonal changes preserve angles and distances in Sym(n). Each point 
in Seal is a fixed point with respect to these changes. Therefore II* is mapped by these 
changes onto itself, the orbit of A resides in II*, and its points are equally-distant from 
A . 

(3) : orbit of each matrix A in M. is an algebraic surface. Because let EV(A) be the 
eigenvalue multi-set for A. By Lemma 1 of Section [2j the orbit for A consists exactly of 
the matrices in M. which eigenvalue multi-set is EV(A). And this latter condition on a 
matrix X is equivalent to the following system of n + 1 equations on X: 

discr(A) = 0, 

Ci(X) = sigO') • elSym i (ev(A)), < % < n - 1. 

These equations, — except the first one, - are the Viete expressions for a polynomial 
coefficients via its roots, in which there are substituted the values for the roots. Here 

C{(X) is the coefficient of degree i in the characteristic polynomial for X - - a 
polynomial in the variables-elements of X, 

ev(A) is the number sequence made of the values of the multi-set EV(A), with 
repetitions of multiple values, 

elSynij is the elementary symmetric polynomial of degree n — i, 

sig(z) is the appropriate sign: plus or minus. 
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Minimal equations algorithm for orbit of a matrix in Ai 
When given the value n and the values a, bi, ... , 6„_2 for the variables A, /ii, . . . , /i„_2, 
a minimal system for the orbit of a matrix in M. having this eigenvalue multi-set is 
obtained as follows. In the algorithm of the Section [3j add the following linear equations 
to the orthogonality conditions on the matrix Y: A = a, H\ = bi, . . . , /i„_2 = & n -2- 
And further, it is applied the same algorithm for finding relations. Evidently, this 
produces a degree-minimal equation system for the orbit. 

Further, the orbit of A is a prime, connected and compact algebraic surface. Because 
the map 

DO A ■ SO(n) -> Orbit (A), DO A (g) = g c A 

is a polynomial map of algebraic surfaces. By definition of an orbit, its image is Orbit (A), 
and the group SO(ri) is compact and connected. Hence, primality of Orbit (A) follows 
from a) that SO(n) is prime, b) that DO a is a polynomial map onto whole Orbit (A), 
c) Lemma A2. 

So: the orbit for A is an algebraic, prime, compact, connected surface contained in 
some sphere inside a hyper-plane II*. The minimal dimension for an embracing plane 
depends on the multiplicities in the spectrum of A. 

The further discourse in this lemma (on a smooth parameterization of an orbit) we 
provide only for the case of a matrix A of a maximal spectrum in M.. This is sufficient 
for the main two theorems in this paper. 

Let Bas = {e\, . . . ,e n } be the basis consisting of the eigenvectors of A, and the 
vectors e\ and €2 belong to the same eigenvalue. Then the stabilizer St (A) is a single- 
parameter subgroup of rotations in the two-dimensional plane Li^e^e^)- It is known 
from the Lie group theory that in this case the orbit of A is locally diffeomorphic to 
a smooth surface — homogeneous space SO' of conjugated classes in SO{n) modulo 
St(-D), having dimension dimSO(n) — 1. 

Due to the action of SO(n), small neighborhoods of all points of the orbit are iso- 
metric. This provides an atlas for a smooth parameterization of the orbit. 

Although a smooth parameterization of the orbit is proved, - - by a reference to 
the general Lie group theory, — still we provide in Application a certain more definite 
description of parameterization for the case of n = 3. 
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5.2 Orbit structure details in the case of n = 3 
Two definitions: 

1) for a point M in an affine space and for a vector subspace L, H(M, L) is the 
plane with the subspace L, containing the point M. 

2) The expression "subspace of 1-orbit" means the vector subspace of the plane of 
1-orbit. 

Lemma ST. 

In the case of n — 3, the following holds for orbit of any non-scalar matrix Am M. 

(1) This orbit resides on the four-dimensional sphere (in a five-dimensional plane) 
and is diffeomorphic to the projective plane. 

This is a bi-dimensional, algebraically-quadratic surface. It resides in the hyper- 
plane orthogonal to the line of scalar matrices and it is not contained in any four- 
dimensional plane. 

This orbit is union of a smooth single-parameter family of circumferences of the same 
radius having a common point A. 

The orbit orthogonally intersects the plane Dg of diagonal matrices, and it has 
exactly three common points with Dg. 

(2) This orbit is described by four quadratic equations in the variables x 2 ,2, ^3,3, 
xi,2, 2:1,3, ^2,3 of the elements of a matrix X in Sym(3). 

(3) The orbit diameter for a matrix in M. having the eigenvalues (A, A, //) is |A — fj,\. 

Proof. Due to the group action and by the statement (7) of Lemma 1, it is sufficient 
to prove Lemma ST for the diagonal matrix D = diag(l, 1, —2). 
Let us prove the statement (1). 

Represent RP(2) as the manifold of axes in the three- dimensional space V, and 
define the map 

s : RP(2) ->■ Orbit(D), s(l) = Rt(e 3 ,l) c D. 

Here 1) D has the eigenvalue —2 along the axis of e^, 

2) I is an arbitrary axis in V, 

3) g = Rt(es, I) is the operator of rotation of V from the line of e% to the line / around 
the normal to the plane L(es, I), and we put that this rotation is identic in the case of 
I = e 3 , 

4) g c is the change in Sym(3) by this rotation. 

This map is a diffeomorphic injection from the projective plane into the space of 
Sym(3), onto the orbit. This holds so due to the following reasons, 
a) each form A in the orbit of D has the eigenvalue multiset {1, 1, —2} and uniquely 
defines the axis / that has in A the eigenvalue —2. b) A defines uniquely the plane 
L 2 orthogonal to /, the operator A has identical restriction on L 2 , the operator g maps 
the plane L(ei,e2) to L 2 . c) Evidently, it holds g c D = A. 

This simple and important consideration about parameterizing this orbit by finding 
the axis for the eigenvalue —2 has been communicated to me by S. Orevkov. 

Note that even though projective plane is not a sub-manifold in V, the above orbit 
- image of the map s — is an algebraic variety in the 5-dimensional space. 



15 



The orbit also has certain remarkable metrical and algebraic properties. Thus, let us 
start with the orbit construction made by the circumference family. 

Consider first the 1-orbit of D by the axis e\. This is a circumference (residing in its 
bi- dimensional plane), which we denote Ciri. This is visible from the explicit formulae 
for a change in D by a rotation around e\. The general matrix of this 1-orbit is 

C1G = efD = 

[[l, o, o ] 

[0, c~2 -2*b~2, -3*c*b ] 
[0, -3*c*b, -2*c~2 +b~2]] , 

where c = cos0, b = sin0. The only diagonal matrices in this circumference 
are D and D 2 = diag(l, — 2, 1). When the angle changes from zero to tt/2, the 
matrix C1G passes (without repetition) the arc from D to D 2 . At this stage, the matrix 
element C1G(2,3) first changes monotonously from zero to —3/2, and then changes 
monotonously back to zero. When the angle changes from tt/2 to n, the matrix C1G 
passes (without repetition) the arc from D 2 to D; this is the "lower" half of the 1-orbit, 
it differs from the upper half only in the sign of the element C1G(2, 3). 

The subspace L\ of the 1-orbit Cirx is linearly generated by the vector D 2 — D and 
the tangent vector to Cirx in the point D. After norming by the factor —3, these two 
vectors produce the following basis for the space L\\ 

Ml = diag(0, 1, -1), M2 = [0, 0, 0] 

[0, 0, 1] 
[0, 1, 0] . 

Now, let / be the axis in the space L(e.\, e 2 ) obtained from e\ by rotation at an angle 
ip around e^, let h be the operator of this rotation. Then the 1-orbit Cir(^) for D by I 
is a circumference containing D and isometric to Cir^ this isometry is the operator h c 
(in the space Sym(3)), and it also maps the plane L\ to the plane L(ip) of the 1-orbit 
Cir(^). This property holds due to that the operator h c conjugates the subgroups G(e 1 ) 
and G(l). Therefore, for these two 1-orbits it holds the equation 

Orbit(Z)D= h c Orbit (e^h^D). 

But in this case, h is a rotation around e 3 , and D = diag(l, 1, — 2). Hence h com- 
mutes with the operator D, and the relation between these two 1-orbits is simplified: 
Orbit (I) D = h c Orbit (e\)D. That is these 1-orbits are conjugated by h. 

So, it appears that the orbit for D contains the union of a smooth single-parameter 
family of circumferences of the same radius, all containing the point D. Let us also 
note that the plane II(^) of 1-orbit in this family turns dependently on ip in various 
directions in a five-dimensional space. We shall see later that Orbit (D) is union of these 
circumferences. 

Let us prove that the orbit of the matrix D is not contained in any four- 
dimensional plane. 

The simplest way to see this is to consider an embracing plane for the orbit of 
D = diag(0,0, 1). Because by the statement (7) of Lemma 1, the latter orbit differs 
metrically from the former only by a homothety by a non-zero factor. 

It is easy to compute that the tangent space to the orbit in the point D is generated 
by the vectors 
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[0 0] [0 1] 

Tl = [0 1] and T2 = [0 0] 
[0 10] [10 0]. 



These are the results of changes in D by infinitely small rotation around e\ and 
around e 2 respectively, they are obtained by commuting the matrix D with two (anti- 
symmetric) matrices from a basis of the Lie algebra so(3). 

Also the orbit contains the matrices D 2 = diag(0, 1,0) and D 3 = diag(l,0,0). 
Respectively, the tangent space in D 3 contains the vector 

[0 l 0] 
b = [l o o] 

[0 0]. 

Denote D 2 l = D 2 -D, D 3 f = D 3 -D. 

If a plane L contains the orbit for D, then L contains the point D and also the 
directions Tl, T2, B, D 2 f — diag(0, 1, — 1), D 3 f — diag(l, 0, — 1). The vectors Tl, 
T2 and B are mutually orthogonal (recall the s-metric). The vector TV is orthogonal 
to Tl, T2 and B for i — 2,3, and the vectors D 2 f and D 3 f constitute a staircase 
matrix of rank two. So: the above five directions are linearly independent and belong 
to each plane containing the Orbit (D). 

Let us prove the Lemma statement about quadratic equations. 

A general algorithm for finding minimal equations for the orbit of D is described 
earlier, in the Section (3) and in the proof for Lemma 2. In the case of n = 3 the 
computation is not expensive. Let us print its result for the matrix 

£> = diag(l,l,-2). 

For any other matrix in AiTL, the orbit geometry differs only in homothety by a non-zero 
factor (by the statement (7) of Lemma 1). 

In this example, the orbit resides in the hyper-plane LI*, defined by the equation 
xll + x22 + x33 = o (and orthogonal to the line Seal), and it resides on the four- 
dimensional sphere having its center in zero. The algorithm derives the following equa- 
tions for this orbit: 

orbitEqs = {xll + x22 + x33, — xll is eliminated 

x23~2 - x22*x33 + x22 + x33 - 1, 
xl3~2 + x22*x33 + x33~2 - x22 - 1, 
xl2*x33 - xl3*x23 - xl2, 
xl2~2 + x22~2 + x22*x33 - x33 - 1} 

- four quadratic equations for the variables x22, x33, xl2, xl3, x23 defining a 
bi-dimensional smooth surface in a five- dimensional plane. 

We need to add: "if the program computes correct" . On the other hand, the degrees 
in this example are small, and we believe, one could verify the computation manually. 

The lemma is proved. 
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5.3 Parametric equations for 1-orbit. 

Analyzing and checking the orbit equations 

It is remarkable that the system for the orbit occurs quadratic. But how do we derive 
from these equations a confirmation for our theorem on the orbit structure, for example, 
the statement about circumferene family? 

Bringing to main axes in some of the above equations and summing some of them 
shows that the orbit is intersection of 1) cylinder (with the element space L(xl2, xl3)) 
over a straight circular cone, 2) cylinder over an ellipsoid, and so on. 

Our impression is that it is difficult to verify the orbit structure in this way. There- 
fore, to verify this computer evaluation result, — and also for the algorithm demonstra- 
tion, — let us apply certain interesting additional checks and computations. In passing, 
we also derive a generic parametric equation system for an 1-orbit. 

Let us find, in what way the equations orbitEqs contain the circumferences having 
a common point D = diag(l, 1, —2) and being the 1-orbits by the axes I — as it is 
described in Lemma ST. 

Define an axis I by a real number k, putting that the vector v of this axis is v — 
el + ke2- So: 1) the axis passes though the right-hand half-plane, 
2) k = corresponds to I = ei, 3) k = oo corresponds to I = e 2 . 

Further, for a negative k, the orbit by the axis l(k) differs from the orbit by l(—k) 
only in the sign of the element X(l,3) of the current matrix. This will be visible from 
the parametric formula for an 1-orbit derived below. 

The orbit equations are derived by applying the algorithm of the Section |3j Further, 
to obtain equations for the 1-orbit by I, it suffices to add the fixed-axis condition to the 
orthogonality conditions on the operator Y in the algorithm. The fixed-axis condition 
is Y(e.\ + ke?) — e\ + ke^ - three linear equations. Here A; is a parameter, and the 
algorithm operates with polynomials having rational functions in k as coefficients. And 
the result equations depend explicitly on the parameter k. 

In this example, the algorithm produces a minimal- degree basis for the ideal of 
relations for the elements X(i,j) of the matrix. And by Lemma ST, we expect to see 
a relation basis defining a circumference in the plane L(l). The algorithm prints the 
following equations for an 1-orbit: 

1-orbitEqs = { xll + x33*k~2/ (k"2+l) + (k~2-l) / (k~2+l) , 
x22 + x33/(k~2+l) + (-k~2+l)/(k~2+l) , 
xl3 + k*x23, 

xl2 - x33*k/(k~2+l) - 2*k/(k~2+l), 

x23~2 + x33~2/(k~2+l) + x33/(k~2+l) - 2/(k~2+l) > 

Hence, the general matrix of 1-orbit is 

[-x33*k~2/(k~2+l) +(l-k~2)/(k~2+l) , x33*k/ (k~2+l) + 2*k/(k~2+l), -k*x23] 

[— -x33/(k~2+l) + (k~2-l)/(k~2+l) , x23 ] 

[— — x33 ] 

Here wild-card denotes the elements under the main diagonal in this symmetric 
matrix. And separating of a complete square in the last equation shows that the variables 
x23 and x33 are bound with the ellipse: 

(k~2+l)*x23~2 + (x33+l/2)-2 = 9/4 (E) . 
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The elements X(i,j) are expressed linearly though x23 and x33. So, we see that 
each 1-orbit is an ellipse which is the projection of the ellipse (E) to the corresponding 
plane. This agrees with the statement of the Lemma ST about the circumference family. 
One could continue this test and make sure that this ellipse is a circumference of the 
same radius square (9/4), as Ciri from the proof of Lemma ST. 

Further, for the equations 1-orbitEqs there are the two remarkable values for k: 
zero and infinity, k — must produce the equations for the circumference Cir!. This 
substitution yields the matrix 

[l, o, o ] 

[0, -x33-l, x23] 

[0, x23 x33] , x23~2 + (x33+l/2)"2 = 9/4 (El). 

This coincides with the parameterization for the circumference Ciri from Lemma 
ST. And the equation (El) must express the equality detX = — 2 = detD. To check 
this, compute det(X) + 2 modulo (El). Indeed, its result is zero. 

With k approaching infinity, I has the limit and there must appear the equations 
on the 1-orbit of D by e2, and this orbit must be the circumference Cir2, of the same 
radius. The equations 1-orbitEqs contain the relation xl3 + k*x23 = 0. For infinite 
k, let us choose xl3 as the independent variable, and express x23 as -xl3/k. With 
this, substituting in the equations 1-orbitEqs of the limit k = +oo produces the 
generic matrix 

[-X33-1, 0, xl3] 
X = [0, 1, ] 

[xl3, 0, x33] , 

with the condition ((k~2+l)/k~2)*xl3~2 + (x33+l/2)~2 = 9/4. This condition is 
equivalent (under the given limit for k) to the equation xl3~2 + (x33+l/2)~2 = 9/4. 
And this coincides with the parameterization for the circumference Cir2 being the 1-orbit 
by the axis e2- 

5.4 Smooth parameterization for the surface AiTi 

Let us prove the statement (5) of Theorem DS. 

Due to the action of the group SO(n), it is sufficient to build a smooth parameter- 
ization of M for a neighborhood of a diagonal matrix D in MH having D(l) = D(2). 
Recall that II = 11(1,2) is the plane of diagonal matrices D in M satisfying the 
condition D(l) = D(2). 

Let us parameterize a neighborhood U of a matrix D in the plane IT in M.D by 
the vector of n eigenvalues of spectrum width n — 1: diag(A, A, /ii, . . . , \i n -2)- Consider 
the polynomial map 

DO : II x SO(n) ->■ MH, 
DO((X, ... , fi n -2),g) = 9° diag(A, A, . . . , //„_ 2 ). 

By the statement (4) of Theorem DS proved earlier, Orbit (D) is a smooth surface of 
dimension dimSO(n) — 1. And its proof includes a local parameterization of this 
orbit by the conjugation classes by the stabilizer St(D), where this stabilizer consists 
of rotations of the bi-dimensional plane L(e±, e^). Consider the map DO g which differs 
from DO in that operator g in the argument is taken not from 5*0(3) but from the 
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space of the conjugation classes by the stabilizer. Fixing the first argument in DO q to 
correspond to D, we obtain a local parameterization for Orbit (D). 

In the neighborhood U the values in the vector (A, /ii, . . . , /i„_ 2 ) remain different. 
Hence, in U all the matrices have the same vector of eigenvalue multiplicities as D: 
(2, 1, ... , 1). Hence, all the matrices in this neighborhood (in n) have the same sta- 
bilizer in SO(n). And it was proved earlier that the orbit intersection with the plane 
of diagonal matrices is orthogonal. Due to all this, it is evident that DO g is a smooth 
parameterization of a neighborhood of the matrix D in M (this is the same as in MM). 

The statement (5) is proved. 

Corollary about dimension of M 

In particular, dimension of M is dim(n) + dim(Orbit(D)) = 

ra-l + dimSO(n) -1 = n+ (n(n - l)/2) - 2 = ra(n + l)/2 - 2 = dim(Sym(n)) - 2. 

Therefore (for each n), co-dimension of .M is 2. 

So, MM is represented as a smooth fibration of co-dimension 2 over a (n — 1)- 
dimensional plane with the fiber being an orbit and a smooth, algebraic, compact, 
connected surface of dimension dim(SO(n)) — 1. 

5.5 Bounds on an orbit diameter 

Let us prove the statement (3.3) of Theorem DS. We need to derive bounds on the 
orbit diameter for a matrix in MM. Such a matrix is diagonalized by some change in 
SO(n), and this change preserves distances. Therefore it is sufficient to prove bounds 
for a diagonal matrix D = diag(Ai, . . . , A n ) in MM. Let % and j be so that the value 
|Aj — A j- 1 is maximal for all the pairs of the eigenvalues of D. Orbit (£)) contains all 
the results of changes in D by rotations of the plane L(e i: e,). With these changes, the 
result matrix passes through a circumference of the radius |Aj — Aj|/2. This proves the 
first inequality of the statement (3.3). 

The upper bound of this statement follows from that (by Lemma 2), the orbit resides 
on the hyper-sphere with the center in the scalar matrix, being the projection of D to the 
scalar line. Hence, the diameter of this sphere is an upper bound on the orbit diameter. 
This proves the statement (3.3). 

A detail: the diameter square of the above sphere is expressed via the eigenvalues 
as (Ai — m) 2 + . . . + (A n — m) 2 , where m is the arithmetical mean value for the numbers 
Ai, • • • , A n . 
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5.6 Proof for the structure of Ai for n = 3 

By statement (8) of Lemma 1, Ai is a straight cylinder over the surface restriction A4q. 
So, it remains to find the geometric structure of M.q. The eigenvalue multi-set of each 
matrix in Aio is of the kind {A, A, — 2A}. By statement (7.2) of Lemma 1, the orbit 
of such a matrix is obtained from the orbit of the matrix D = diag(l, 1, —2) by the 
homothety by the factor A. Therefore, Aio is the cone with its vertex in zero and its 
base being the orbit Orbit (D). This (bi-dimensional) diffeomorphic image of projective 
plane (Section Lemma ST) resides on a four-dimensional sphere 5*4 with center in 
zero. 

Any diagonal matrix belonging to Aio is of one of the following kinds: -D(A) = 
diag(A, A, — 2A) - and also the falilies D 2 (X) and -D 3 (A) obtained from D(X) by per- 
mutations on the main diagonal. 

Remarks: 1) For each A, the matrices D(X), D2(X) and D%(X) belong to the 
same orbit; 2) D(l) and D(—l) are central-symmetric points on the sphere S4, but 
they belong to different orbits; 3) Intersection of Alo with the sphere S4 is union of 
Orbit (D(l) and Orbit(D(-l)). 

Theorem. 

Zero matrix is the only singular point on the surface Aio- 

Proof. Regularity of any non-zero matrix in Aio is proved earlier, in Subsection 

El 

Why the vertex is singular? As the orbit of D surrounds the vertex in a certain 
complex way, we cannot simply see singularity of the vertex by referring to the picture 
of a straight circular cone. Let us prove singularity by considering linear dependency 

between the tangent vectors. 

The vertex (0) is singular because 1) Aio is a cone with vertex in zero, 
2) we present four (straight) lines, intersecting in zero, belonging to Aio, and having 
linearly independent directions. The two of them are the lines, connecting the point 
with the matrices D = diag(l, 1, —2) and D 2 = diag(l, —2, 1) respectively. Further, the 
orbit of the matrix D contains the 1-orbit for D by the axis e\. This is a circumference 
containing the matrices D and D 2 , so that D corresponds to the angle 0, and D 2 
corresponds the angle tt/2, and these two are the only diagonal matrices in this 1-orbit. 

For the proof, we also need some non-diagonal matrices. For example, e[^ A D is the 
matrix 

[10 0] 
Ml = [0 a b ] 

[0 b -a-1] , a = -1/2, b = -3/2. 

Similarly, the 1-orbit for D by e 2 contains the matrix 

[a b ] 
M2 = [010] 
[b -a-1] 

- with the same a and b. It is easier to see linear independence of the sym- 
metric matrices D, D 2 , Ml, M2 if we represent each of them in the row-vector form: 
[X(i, j)| 1 < % < j < 3] - with skipping the sub-diagonal part of the matrix. This 
produces the matrix of the four rows: 
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[1, 0, 0, 1, 0, -2 ] 

[1, 0, 0, 0, -2, 1 ] 

[1, 0, 0, a, b, -a-1] 

[a, 0, b, 1, 0, -a-1]. 

The first loop of clearing of the first column in the staircase form algorithm produces 
the matrix 

[1. 0, 0, 1, 0, -2 ] 
[0, 0, 0, -1, -2, 3 ] 
[0, 0, 0, a-1, b, -a+1] 
[0, 0, b, 1-a, 0, a-1 ] 

Further, the fourth row is moved to the second place, and to the third row it is added 
the second row multiplied by a-1. This produces the matrix 

[l, o, o, l, o, -2 ] 

[0, 0, b, 1-a, 0, a-1 ] 
[0, 0, 0, -1, -2, 3 ] 
[0 , , , , b2 a2 ] . 

Here b ^ 0, b2 = b - 2 (a-1) = -3/2 - 2 (-3/2) ^ 0. Hence, the matrices 
D, D2, Ml, M2 are not linearly dependent. So, we have found such four straight lines 
in the surface A^o which are linearly independent and intersect in zero. Therefore, zero 
is a singular point in A4q. Let us explain why it is singular. All other points of this 
surface are regular, and the tangent space in each of non-zero points has dimension 
three. To prove by contradiction, suppose that zero in regular. Then, by the definition 
of a regular point, some neighborhood of zero in Ai is diffeomorphic to a ball in M 3 . 
Such a ball cannot contain four linearly independent directions which we have found. 
This contradiction shows that zero is a singular point in .Mo- 
Theorem DS3 is proved. 
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6 Application 



6.1 Several lemmata 

Lemma AO (probably, known). 

For a diagonal quadratic matrix D and an anti-symmetric matrix A of the same size, 
the commutator [D, A] is a symmetric matrix having zero main diagonal. 

Proof. Symmetry of the result matrix is proved by the following equations: 

(DA - AD)* = (DA)* - (AD)* = A*D* - D*A* = (-A)D - D(-A) = DA- AD. 

Further, denote Di the diagonal matrix in which D(j) = 5 it j for all j. Each diagonal 
matrix is a linear combination of the matrices D^. And as commutator is a bi-linear 
map, it is sufficient to prove the goal for the commutator [A, A] for each i. The given 
matrix A has zero main diagonal, and it is visible that A ■ D{ and Di ■ A have zero 
main diagonal. Therefore, their sum has zero main diagonal. 
Example: 

10 1 1 * | a 01 - |a 01*10 1| = |0 b|-|0 a I = |0 b-a| 
1-1 01 10 b| 10 b| 1-1 01 l-a 01 I -b 1 Ib-aO I. 

Lemma Al (known). 

For any polynomial map F : R n — > R m , counter-image Ai' of an algebraic set Ai 
is an algebraic set. 

Proof. Let Ai be the set of zeroes of polynomials {p±, . . . ,Pk}- Then Ai' is the set 
of zeroes for the set of compositions {p 1 (F(X)), . . . ,p^.(F(X))}. Indeed, if A belongs 
to Ai', then F(A) belongs to Ai; then pi(F(A)) = for all i. 

Conversely: if pi(F(A)) = for all i, then F(A) belongs to Ai, hence A belongs to 
Ai'. 

Lemma A2 (known). 

For a surjective polynomial map from an algebraic set Ai onto an algebraic set Ai', 
if Ai is prime, then Ai' is prime. 

To prove by contradiction, suppose that Ai' is union of algebraic sets Mi and M2, 
none of which contains another. By Lemma Al, the counter-images T\ and T2 of M\ 
and of M 2 respectively are algebraic sets. Also it holds Ai = T x U T 2 . And as none of 
Mi and M 2 contains another, their counter-images are in the same relation. This proves 
the lemma. 

6.2 Addition to Lemma 2. 

a) St (A) is not a normal subgroup, and a linear complement in the Lie algebra to the 
tangent space for this subgroup is not closed by commutator. 

b) Although a smooth parameterization is already proved (by a general reference to 
the Lie group theory), still let us describe a local smooth parameterization by a more 
definite construct, and for simplicity, consider only the case of n = 3. 

Due to translations by 5*0(3) on the orbit, it is sufficient to define a smooth pa- 
rameterization for an orbit of a diagonal matrix D = diag(a, a,b), with a ^ b, in a 
neighborhood of D. 
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In this case, the stabilizer St(-D) is the subgroup of rotations around e^. Let X3 be the 
tangent vector in unity to St(D), let Xi an d X2 be the tangent vectors in unity to the 
(uni-parametric) subgroups of rotations around e\ and e 2 respectively. The sub-space 
V = L(xi, X2) does not contain the vector X3 (but it is not closed by commutator). 

Let expy be the map exp of the operator exponent considered in restriction to 
V. expy maps V to a bi-dimensional surface inside SO(3). There is also the map 
C : 5*0(3) — ?■ Sym(3) of the group action: C(g) = g c D. 

The composition CE(A) = C(exp v (A)) is a map from V to Sym(3), having image 
inside M.. It also is a composition of a pair of smooth maps and has rank two in zero. 
The latter holds due to the following two reasons. 

1) The infinitely small rotation around e\ affects only the rows No 2 and 3 in the matrix 
D, and in the tangent map, the image element at the position (1,3) is zero, and the 
element at (2, 3) is b — a. 

2) The infinitely small rotation around e 2 affects only the rows No 1 and 3 in the matrix 
D, and in the tangent map, the image element at the position (1, 3) is b — a, and the 
element at (2, 3) is zero. 

So, this produces a diffeomorphism from a neighborhood V e in a bi-dimensional plane 
onto a neighborhood in the orbit of D. 

For the case of arbitrary n, — in a similar way, — the tangent vector to the stabilizer 
has a linear complement V of dimension dimSO(n) — 1, and there is considered the 
parameterization map A — > T(exp v (A)). 
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